Mixing of axino and goldstino, and axino mass 
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Axino, related to the SUSY transformation of axion, can mix with goldstino in principle. This 
case is realized when some superfields carrying nonvanishing Peccei-Quinn charges develop both 
scalar VEVs and F-terms. In this case, we present a proper definition of axion and axino. With 
this definition, we present the QCD axino mass in the most general framework, including non- 
minimal Kahler potential. The axino mass is known to have a hierarchical mass structure depending 
on accidental symmetries. With only one axino, if Ga = where G = K + \n\W\ 2 , we obtain 
ma = m 3 / 2 - For Ga 0, the axino mass depends on the details of the Kahler potential. In the 
gauge mediation scenario, the gaugino mass is the dominant axino mass parameter. Therefore, we 
can take the theoretical QCD axino mass as a free parameter in the study of its cosmological effects, 
ranging from eV to multi-TeV scales, without a present knowledge on its ultraviolet completion. 
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I. INTRODUCTION 

The lightest supersymmetric particle (LSP) with Im- 
parity conservation is absolutely stable and can con- 
tribute to the present energy density of the universe as a 
dominant component of cold dark matter (CDM). If the 
QCD axion solves the strong CP problem, its fcrmionic 
SUSY partner, axino, must be considered in the CDM 
estimate in cosmology. Firstly it can be a natural candi- 
date for CDM if it is the LSP [1,11]. Second, if the axino is 
much heavier than the neutralino LSP, still it affects the 
mass density estimate of the neutralino LSP because the 
nonthermal decay of axino may dominate the estimate 
Therefore, it is of utmost importance to clarify what 
is the axino mass in a specific supcrgravity model. The 
axino mass is obtained when SUSY is broken. A naive 
guess on the axino mass is of order m 3 / 2 since the axino 
mass is the soft term. But a leading loop corrections can 
reduce it to 01113/2. If this is suppressed, the next level 
hierarchical masses are arising from the gauge mediation 
and even some accidental symmetries can reduce it fur- 
ther to m s (fa/M P ), m s (/ 2 /Mj,), ■■■ @. In fact, in the 
literature, axino mass has been considered in a vast mass 
regions from eV to trans- TeV [Tl— f]~5| . 

The axino mass depends on two symmetry breakings, 
the Peccei-Quinn (PQ) symmetry breaking and SUSY 
breaking. The effect of the PQ symmetry breaking must 
introduce a massless axion except from the contribution 
of the anomaly term and the effect of SUSY breaking 
introduces a mass parameter m 3 / 2 . The contribution to 
the axino mass from the parameter 7713/2 is arising by 
the F-tcrms while the PQ symmetry breaking is given by 
the VEVs of scalar fields. In this paper, we study these 
two contributions in a most general form, and express 
the axino mass in terms of 7713/2 with the general Kahler 
form. 

In global SUSY models, the axion a is defined through 
the Peccei-Quinn(PQ) transformation for Vdw = 1 fields 



fa6 ia/f a =Y J v i e u 

i 

Ja 



(1) 



where Ti,Vi and at are the eigenvalue of the PQ charge 
operator T, the vacuum expectation value(VEV) and the 
phase of <f>i, respectively. The PQ direction of <Zj is TiO = 

A prototype axion model with a global SUSY needs 
at least three chiral fields, to introduce a VEV breaking 
the PQ symmetry fl6l | . In this introduction, we do not 
introduce SUSY breaking, but only introduce the PQ 
symmetry and its breaking. The supcrpotcntial having a 
global PQ symmetry is 



W = R(S 1 S 2 -f 2 a ), 



(2) 



where the PQ charges of R, Si, and S2 are 0, +1, and — 1, 
respectively. Equation ([2]) has an additional i?-symmetry 
whose charges are 2, 0, and 0, respectively for R,Si, 
and 52- To have the standard kinetic energy terms, the 
Kahler potential is taken as K = RR* + SiSf + S 2 S|. 
Then, the potential V is given by 

V = \S 1 S 2 -f 2 \ 2 + \RS 2 \ 2 + \RS 1 \ 2 . (3) 

which is minimized at (S1S2) = f 2 and (R) = 0. To 
show the superTrace (STr) of M 2 , we choose the fields 
near (Si) = (S 2 ) = f a , 



51 = -j=(\/2f a + pi + ia-y) 

5 2 = -j=(y/2f a + P2+ ia 2 ) 

R = ^(/°r + ia R)> 



(4) 
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where p's are scalars and a's are pscudoscalars. The 
mass matrix for CP even scalars (ppt, pi, P2) and CP odd 
scalars (or, 01,02) are given by the same squared mass 
matrix, 



(5) 



so eigenvalues are given by (2/ 2 , 2/ 2 , 0). One CP odd 
scalar should be massless since it is the goldstone boson 
corresponding to the spontaneously broken PQ symme- 
try. 

On the other hand, fermion mass matrix is given by 



fermiomc partner axino a, 





(6) 



whose eigenvalues are (y/2f a , —\2f a , 0). Supersymmetry 
is not broken since (V) = 0, so should produce STr M 2 = 
0, 



STrM 2 



'CPevcn 



'CPodd 



2m 



fcrmion 



= 0. 



(7) 



In Eq. ([5]), the smaller mass is the axino mass, and one 
of the larger ones is the hypothetical goldstino mass. This 
larger one becomes exactly massless when the supcrHiggs 
mechanism is operative as will be discussed later. 

The superpotential Eq. @ can be rewritten in terms 
of R, the superfields S = (Si + S^Vv^, and the axion su- 
pcrficld A = (Si —S2)/ v2 containing axion (01 — 02)/ v2j 



W = R 



S 2 -A 2 



r 



0. 



(8) 



The VEVs are (R) = (A) = and (S) = y/2f a . After 
integrating out R, there does not exist low energy self 
interactions of A due to the nonrcnormalization theorem. 

Here note that the zero mass eigenstates of scalar, 
pscudoscalar and fermion indicate that they are related 
by SUSY transformation and form the axion supermul- 
tiplet A. In general, the same mass eigenstates of scalar, 
pscudoscalar and fermion are not related by supersym- 
metry transformation and do not form a supermultiplet. 
For the supermultiplet condition, interactions must be 
supersymmctric also. In our case, however, the zero mass 
eigenstates form a supermultiplet A which survives down 
to the low energy scale. This renders the nonlinear rep- 
resentations for Si = Lpe A l$ a and S2 = ipe~ A /f a , which 
show explicitly the shift symmetry of the axion super- 
field. 



II. THE PQ SYMMETRY IN SUPERGRAVITY 

The supersymmetrization of axion models introduces 
a full axion supermultiplet A which contains the pseu- 
doscalar axion a, its scalar partner saxion s, and their 



1 



A=—(s 



la ) + s/2M + F A &&, 



(9) 



where Fa stands for an auxiliary field and •& for a Grass- 
mann coordinate. 

In the supersymmctric version of axionic models, the 
interactions of the saxion and the axino with matter are 
related by supersymmetry to those of the axion. Sax- 
ion and axino are better to accompany a in Eq. (JlJ 
to preserve SUSY. In Eq. (pQ), f a is a VEV of some 
real scalar field which is called pp. Supersymmetriza- 
tion of ifn needs its pseudoscalar partner tpj, forming a 
complex scalar (p. When the PQ symmetry is not bro- 
ken, the PQ charged fields are of the ip type, its phase 
changes when the PQ transformation is performed, and 
the fields are the real and the imaginary components of 
(p. On the other hand, if it is spontaneously broken, the 
PQ charge is not realized unitarilly but realized in the 
Nambu- Goldstone manner, i.e. the Goldstone boson a is 
created, goes up to the phase as in Eq. ([1]), and the PQ 
symmetry is its shift symmetry. [The PQ symmetry is 
nonlinearly realized on the action to a.] A nonzero SUSY 
breaking F-term of the <p type fields signals the PQ sym- 
metry breaking also if it carry PQ charges. Therefore, the 
direction of axion does not necessarily coincide with the 
direction of axino which is going to be orthogonal to the 
goldstino which in turn is determined by the F-tcrms. 
In this paper, we define the axion and axino properly, 
and set up the formulae for the axino mass even in case 
that a mixing of axino with goldstino is present. Here, 
the axion component is still defined by the coefficient of 
■& term since the F-term or the coefficient of $ 2 term is 
auxilliary. 



A. Origin of axino-goldstino mixing 

The PQ symmetry (as any global symmetries) in su- 
persgravity has a meaning if both the superpotential W 
and the Kahler potential K respect the symmetry. Ex- 
pansion of the Kahler potential in powers of 1 /Mp leads 
to the following type, 



K = J2fi({M)gj({^}) + h.c. 



1, j 



(10) 



where {4>i} and {</>J} are sets of holomorphic and anti- 
holomorphic fields, respectively. From the Kahler poten- 
tial, one can obtain its contribution to V as 



V G 



f d 2 d j d 2 dJ2fi(m)gA{^}) + ^- (11) 

J J I.J 



where and are the superfields corresponding 

to {<j>i} and {<^}, respectively. 

Consider the leading term beyond the minimal term in 
K, for example (l/Mp)H u H^X* . This must preserve the 
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PQ symmetry so that X carries the PQ charge Y(X) as 
the sum of T(H U ) and r(Hd). Namely, the PQ symmetry 
is also broken by an F term, i.e. X F of X* This 
can be obtained from the superpotential as 



V 



W ~ XiX 2 X 



(12) 



where T(X 1 )+T{X 2 ) = -(T(H U ) + T(H J t)). In this case, 
ji = —X\X 2 /Mp is obtained in Ref. [20(. This example 
shows that it is sufficient to consider the PQ charge car- 
rying scalar components to pick up the axion component, 
and should not include the PQ charge carrying F-terms 
for a definition of the axion. Otherwise, we double count 
some components. 
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FIG. 1: The definition of 7V D w 
phase. Here, JVdw = 4. 



in the Nambu-Goldstone 



1. Introduction of tp type fields 

In non-supcrsymmctric case, as axion is defined in Eq. 
dU, and its property determined by the U(1)pq symme- 
try. In the Wigner-Weyl (WW) realization of the PQ 
symmetry, the PQ charged fields transform as 



(13) 



where I\ = -id/86. In the Nambu-Goldstone (NG) 
phase, there appears a Goldstone boson a which can be 
a combination of the original phase fields. Let a be in $ 
such that in the WW phase it is expanded as 



<I> 



y^^cjXj 



(14) 



In the NG phase, the probability amplitude for a to be 
in the phase of Xi is Cj = Vi/V a where V a = ($) and 
Vi = (Xi). The PQ operation on Eq. flU) is 



(15) 



When the PQ symmetry is realized in the NG man- 
ner by giving Xi its VEV Vi, then the charge operator 
r is not unitarilly realized, and then we must use the 
shift symmetry of the phase fields a, instead of the PQ 
symmetry operation T, with the original information on 
the eigenvalues Ti. So, we apply the infinitesimal shift 
symmetry on ^ qX;, 5a = f a o T 9 and Sai = Tif a 59, and 
we use the relation proportional to 89 to code the spon- 
taneous symmetry breaking. For $ ~ e * r W/<< wnere j a 
is determined by the axion-gluon-gluon anomaly term, 
then by acting the shift operation on $ of Eq. (|15l) in 
the NG phase we obtain r a $ as ^Avi/VajTiXi where 
c» is i>i/V a which is the probability amplitude for a to 
be in X t . In 5a t = Tif a S9, Ti encodes the original de- 
generacy number of the phase of Xi as the axion field 
a completes one period in the NG phase. Representing 
* = (V a +p±)e ir " a/fa and!, = (v i +p i ±)e iria / Ja , where 
p± is in the perpendicular direction to axion's Mexican 



hat valley and has mass proportional to V a ■ Since the 
Pi± masses are different, they should belong to different 
supcrficlds. From this discussion, we obtain the axion 
superfield A in 



A/fa 



(16) 



where Cj = Vi/V a is used. Note that we used e A ^ a on 
both sides in view of Fig. Q] since the axion shift of 27r/ a 
is fully accounted for by e A ^ a . Any integer cannot be 
multiplied or divided in the exponent. Note that ipA 
is composed of two real fields p± and Im tpA , so is ifi , 
and their VEVs are (<pa) = V a and (ipi) = Wj. Also, A 
is composed of two real fields s and a, and its VEV is 
vanishing (A) = 0. 

In the NG phase, we must state the domain wall num- 
ber by the axion shift. In Fig. [T] A^dw = 4 is schemati- 
cally shown. The decay constant f a is given by the coef- 
ficient of axion-gluon-gluon anomaly [lTj. So, the axion 
potential returns to itself by a shift of 2nf a [l8| . The 
original domain wall number is given in Fig. [T] by the 
length (ip). 

Functions of the Nambu-Goldstone fields are multival- 
ued. In the original field space, the field returns to itself 
after the a shift of 2irNjy^j f a . If two axion directions 
have two domain wall numbers n\ and n 2 , the multiplic- 
ity of the vacua is the least common divisor of n\ and n 2 
[181 ]. Thus, we can write ri\ = iVowA and n 2 = Adw/^2, 
where (3i and /3 2 are relatively prime. So, Eq. (|16j) can 
be written as 



1 



Adw ^ A 
Therefore, we obtain 



A/fa = 



E 



i 



V a AW/3 



■<fi e 



A/fa 



HiUilPi 



(17) 



(18) 



whose literary form has the domain wall number 1 since 
do not have a common divisor. From (|17j) . we obtain 



v 

' n 



(19) 
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Now, consider the special case Eq. ([2]) where in 
terms of three fields an ultraviolet completion is achieved, 
W = R{SiS 2 - fa) [II- With the complete knowledge 
on W of Eq. @, A and tp are obtained in terms of Si 



and S 2 as Si = pe A ^ a and S 2 = tpe 



-A/f a 



Writing 



the low energy field A in terms of the high energy fields 
R, Si and 52 is not of much use in the region where only 
the low energy effective fields and the PQ quantum num- 
bers of the original fields are known. In most cases, the 
information on tp type fields is not needed. 



2. Appearance of A type fields m W 

For the Kahler potential, we can have any function of 
A + A* . But, a supersymmetric W with A is 



W(A) = 0. 



(20) 



The U(1)pq invariance guarantees that the axion super- 
field A does not appear in the superpotential. But, below 
the SUSY breaking scale we can introduce the soft terms 
in W by introducing the auxilliary field and respecting 
the shift symmetry, 



W = M 3 Qe aA/f ° 
9 = 1 + msusY"^ 2 - 



(21) 



where M is a parameter and TOsusy = fis is the pa- 
rameter describing the SUSY breaking soft terms. Then, 
ln|Fy| 2 appearing in local asusy is invariant under the 
shift of a, with A defined in Eq. (jTBJ). 

If there are more spontaneously broken global 
U(l) J 4 i (i = 1,2, •••) symmetries, coaxions ai with the 
decay constant /j must respect the shift symmetries and 
we must consider the following W, 



W = M 3 Oe aA/f « 



iAi/fi 



the one orthogonal to it becomes heavy. So, at low en- 
ergy, we have the exponential field A and the tp type field 
X ew . X ew does not accompany a phase field. The coef- 
ficient of exponential of A is very large, and A does not 
appear explicitly in the superpotential. Axino and sax- 
ion are in A. So after integrating out our double counted 
fields, we end up with a supcrficld X cw and a supcrficld 
A. So W of Eq. @ has the same degrees at low energy 
in the NG phase. This is the way to write down the low 
energy theory corresponding to Eq. fl2J). How X cw cou- 
ples to the other light fields depends on the ultraviolet 
completion. Note also that the axion interaction depends 
on axion models [25|, . 

On the other hand, Ref. [24[ considered tp,tp 2 : tp 3 and 
tpH u Hd terms without the \x term. 



4- Comments on the model-independent axion 

The model-independent axion in superstring models is 
combined with the dilaton to make a supermultiplet 27j , 



D = \ 

g 2 



+ i 



ttMI , , iMlJaui/fMi 

8irM P ' 8tt 



(23) 



(22) 



where f M i ~ 10 16 GeV [28|], and (s) ~ 2Mp is not the 
tp type field. Because the corresponding U(l) is gauged, 
qmi is absorbed to the U(l) gauge boson, and the U(l) 
symmetry remains as a global PQ symmetry below the 
scale /mj- Only for this anomalous model-independent 
axion in string models, there is no accompanying tp type 
field. Below fui, the resulting pseudo-Goldstone boson 
will accompany a tp type field. We speculate that this 
model independent axion is the only place for the axion 
not accompanying its tp type field. 



B. Goldstino, axion and axino 



3. Comments on NpqMSSM 

In the literature, models based on NpqMSSM have 
been considered 0, [HJ . 

In Ref. [HI, the original fields were used to show the 
existence of light field X ew which must correspond to 
our tp type field in the effective low energy theory frame- 
work. Also, by considering the original fields in a specific 
model, it was shown that tp, tp 2 and tp 3 are not present. 
In this effective theory, the PQ symmetry is broken, and 
the tp type fields do not carry the PQ quantum number. 
So, (/i + tp)H u Hd are the allowed interaction. At the 
fundamental level, the /x term arises from the ori gina l 
fundamental fields, dictated by the PQ symmetry |20| . 
Counting the number of degrees of freedom, we double 
the fields Si and S2 to tpie Al and tp 2 e A2 . As commented 
before, the light A type field is A = A x - A 2 . Ai+A 2 be- 
comes heavy. For the tp type fields, only X ew is light and 



The axion component is defined in Eq. (fTJ). So, what- 
ever the non-vanishing PQ charge carrying F-terms are, 
the axion is properly defined only by the PQ charge carry- 
ing 1? terms. However, the nonvanishing F-terms define 
the goldstino component. 

Supersymmetry is spontaneously broken when the po- 
tential has nonzero VEV, (V) = J2i > where 
F l = K^F-y Then, there should be a masslcss fermion, 
goldstino. In supergravity, it is absorbed to the longitu- 
dinal component of gravitino ipa through the super- Higgs 
mechanism. The goldstino supcrficld, to which goldstino 
belongs, can be defined by 



F 



(24) 



where F = F l Fi which becomes the F-term of Z . 

Among Xi, the axion superfield is defined by the PQ 
charges of Xi. All the other chiral fields orthogonal to A 
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F 1 F 2 F 3 



t 



X2 



t 

spin-O: \X]\ 



A", 




FIG. 2: The case for more than one axino. It is ambiguous to 
choose the partner of the QCD axion. Normalization toward 
the canonical kinetic term is not depicted in the figure. 



are called coaxino directions. Then, we can consider two 
cases in which the axion supcrficld A allows 

• F A ^ 0, or 

• Fa = 0, but F A ^ from Kahler mixing with other 
SUSY breaking fields. 

This case is shown in Fig. [2j 

With the canonical Kahler potential, F l is just F* , and 
the goldstino fermion gives the SUSY breaking direction 
exactly. On the other hand, in the presence of a Kahler 
mixing, F % is nonzero even though Fi = and seems 
to contribute to the goldstino. However, the physical 
goldstino, defined by the zero mass eigenstate with the 
canonical kinetic term, does not contain such a state. 
This can be shown as follows. In supergravity, fermion 
kinetic and mass terms are given by 



(25) 



where is a general covariant derivatives and m'^ = 
m 3/2[ViGj + (l/3Mp)GiGj] is the fermion mass with 
the goldstino field moded out. The primed fields and 
the primed mass matrix are in the interaction basis. To 
obtain the physical states, we first make the kinetic terms 
canonical and then diagonalize the mass matrix. With 
the redefinition of ip n = V 1 aip a , Eq. (|25p is written as 

+ \[V T m'V] ab r^ b + \[v ] m ,] v*uri> 1 + ■■■ 



l.a.i.b i If .i.a.i.b 



(26) 



Requiring the canonical kinetic terms as V^a^G^V 1 



S a b, or V T b l G il V* 



Sbs,, the mass in this basis is 



m ab = [U T m'U] Q b which we have to diagonalize. 

The matrix V can be written in the form of V — US 
where U is a unitary matrix and S is a scaling matrix. 



Taking the inverse of the canonical kinetic term condi- 
tion, V T b l G i jV* j a — 5ba, we obtain 



r* — l\a 



(V*- 1 ) n ] G : > i (V T - 1 ) b 
= (S^ljTfjG^i^S- 1 )^ = 5 ab . 

Taking the complex conjugation, we obtain 



(S- 1 U*) a j G* i (US- 1 } i i = 5 al 



or equivalently, 



{s- 1 u^) a jG fi = (stfy 



(27) 



(28) 



(29) 



On the other hand, has zero eigenvalue in 

the direction of G % = G^Gj. Since m' = 
(V T )- 1 mV- 1 = (C/*S*- 1 )?7i(5'- 1 C/t), m^Cp = implies 

m ab (S- 1 W) b jGJ = 0. So, (S^W f jG^ k G- k is the gold- 
stino direction. Plugging Eq. (|2"9"|) into this, we obtain 



(S^l^f jG fk G- k = {SU^Gj 



(30) 



We have rotated fermion as ip' 1 = V 1 a ip a to make the 
kinetic term canonical. So, in the new (physical) basis, 
= {V~ 1 ) a = {S- 1 U^) a i{i>')\ Note that S is a 

real scaling diagonal matrix. Rotation of the direction is 
determined by U. SUSY breaking direction Gi defined in 
the interation basis ip' is rotated by U in the new basis ip. 

Then U^ a 1 G\ is just supersymmetry breaking direction of 
G7 h not G\ 

As indicated within the yellow box in Fig. [21 there 
is an ambiguity in identifying the mass eigenstate axino 
corresponding to the QCD axion. 

When there exists only one axino, there is no ambiguity 
in identifying the mass eigenstate QCD axino, because it 
must be the orthogonal state to the goldstino. For two 
fermions, we can consider two cases separately as shown 
in Fig. E] for Fa = and Fig. @] for F A ^ 0. Even 
for one axino and goldstino case, the one beyond the 
axion multiplct is called the coaxino as indicated in the 
left-hand side of Figs. [3] and [4] with a thick brown bar. 
These relatively simple cases will be discussed explicitly 
in Sec. IIII Al where we calculate the axino mass by the 
gravity mediation. In the regrouping of Figs. [3] and 31 
we violated the supcrmultiplct condition, QA = a, and 
may make an error in estimating the axino mass by an 
order 0(m^ 2 /Mp). 

The goldstino multiplet Z is defined by the total F* 
term, 



F A t-^ F l F A F c 
Z = A + V — X l = A + C, 



i=£A 



(31) 



where C is the sum of SUSY breaking coaxino, orthogo- 
nal to the axino supcrficld. Their bosonic and fcrmionic 
components are not the components of the original su- 
pcrficlds. Since axion and goldstino arc finally defined 
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\F A ^0, Fcyo] -»| F= x/ET^ 



Fa = 



t 
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\s + za 
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F A , ji 
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s' + ia' 



A' 



G(goldstino) 



Z 



\F A ^0, F c ^0\—>- F= v'E, F'Fi 



\Fa ^ 0| L5bJ 

t t 
a'(axino') , c(coaxino) 

t t 
I s + ia I c 



F 4 , ^ 



1 

a 


_L 





G(goldstino) 



s + ia' g 



A' 



FIG. 3: The case Fa = for the axion (blue) and goldstino 
(red) multiplets. The axion direction a is defined by the PQ 
symmetry through Eq. Q and the goldstino (G) and axino 
(a) directions are defined by the fermion mass eigenvalues. 
Normalization toward the canonical kinetic term is not de- 
picted in the figure. 



after SUSY is broken, the mismatch is generic. In view 
of Figs. [3] and U introducing one coaxino, we can regroup 
the axion and coaxino multiplets as 



A=(A,a',F A ) 
C=(c',c,F c ) 



(32) 



where the scalar component of the axion multiplct is A = 
s' + ia. a' in Figs. [3] (and also in U) is not a mass 
eigenstate since the identification of a _L G does not care 
about its supcrpartncr a'. The scalar potential is given 
by 



V 



Mle G ' M i 



[G^GiGj - 3M|], 



(33) 



and pseudoscalar a in A does not appear in V since G is a 
function of (1/2)(A+A*) = s'. It is a massless eigenstate, 
and therefore the pseudoscalar mass matrix is already di- 
agonalized in the left-hand side of the brown bar. On the 
other hand, the fermion mass matrix should be diagonal- 
ized in the right-hand side of the brown bar. One eigen- 
value in the direction of F l , (F A ,F C ) is massless along 
the direction of Z which is goldstino. After making the 
fermion kinetic term canonical, goldstino indicates the 
direction (Fa,Fc). The remaining eigenvalue is inter- 
preted as the axino mass. If Fa <C Fq, 5~ a' — (F A /F)Z 
and the mass eigenstate is axino-like. 

In general, however, after the SUSY breaking, the mass 
eigenstates are 



Scalar 
Pseudoscalar 
Fermion 



s, Re c 
a, Im c 

a = a' — 



= Im c' (34) 
(F A /F)Z, Z, 



where s and c are the mass eigenstates after diagonalizing 
the mass matrix in the (s', c') basis. In the last step, the 
shift symmetry of a must be invoked to guarantee the 
massless axion. 



FIG. 4: The same as Fig. [3j except for F A ^ 0. 



This process makes sense when we consider the axion 
interaction with gluino 



1 

fa 



d 2 $AW a W a 



(35) 



If Fa = as in Fig. [3l the gluino mass does not get a con- 
tribution from the axion multiplet. However, if Fa 0, 
the gaugino mass should be studied carefully in a spe- 
cific model. This is true also if the Kahlcr potential has 
a large mixing term. 

In our language, after the goldstino G is absorbed to 
ijjp, we can write the ?/v interaction as[2l|. 



2V3Mj 



(36) 



So, before hiding the goldstino into ip^, any chiral field 
interaction in supergravity has a coupling suppressed by 
I/Mp, 



(coupling) $ ■ 



(37) 



Therefore, the coupling of A 1 is suppressed by Mp, not 
by f a , in Eq. (|35|) . After absorbing the goldstino into 
ipp, any other chiral fields are orthogonalized not to have 
an F-term but its supergravity coupling is suppressed by 
1/Mp again, and the goldstino supcrficld G couples to a 
gauge supermultiplet as 



Mp J 



(38) 



where Fq = F gives the gaugino mass. In addition, the 
axino interaction with gluino is 



1 

fa 



d 2 $AW a W a 



1 



->• T (2- e a G) (gluino) a (gluon) a 

Ja 

-> -j- a (gluino) a (gluon) a 

Ja 

(^Vv)(gluino) a (gluon) Q 



(39) 



faMp 



7 



a 




Q % 



a 



y 

(S) 

FIG. 5: The two-loop axino mass in the KSVZ model. 



where m 1 / 2 . a are the gaugino masses. These contribu- 
tions give masses of order up to probably 10 GeV. 

In the gravity mediation scenario, there exists another 
parameter m 3 / 2 which can be much larger than 10 GeV. 
Without the axino-goldstino mixing in the Kahler poten- 
tial, the axino direction is the same as that of axion and 
the superpotential determines the axino mass. It means 
that axino mass arises from loop diagrams as in the Fig. 
[5] example. Therefore, without the axino-goldstino mix- 
ing in the Kahler potential, axino mass is not going to 
be larger than 10 GeV. Thus, a very heavy axino mass is 
possible only if there is a significant A — Z mixing in the 
Kahler potential. 

If we consider the anomaly mediation, there would 
be additional contributions from breaking the conformal 
symmetry which appear as anomalous dimensions in the 
superpotential terms [29|. Generally, they are propor- 
tional to 7anom 1^3/2 for the term Wi in the superpo- 
tential. Thus, the axino mass has contributions from all 
the above cases and expressed as 



So axino interaction implies some accompanying grav- 
itino interaction, suppressed by the product f a Mp. This 
makes sense since it must respect the spontaneous sym- 
metry breaking suppressed by f a and the supcr-Higgs 
mechanism suppressed by Mp. 



m a — I ^goldstino + Cf° m J 

\ /—terms in W / 

+ &i TO l/2,a- 



(43) 



a— gaugino 



C. Axino mass with SUSY diagrams 



We expect that £| nom is of order O(10 2 ). An example 
has been discussed in 1301. 



In Eq. ©, the smaller mass was shown as the ax- 
ino mass, and one of the larger ones is the hypothetical 
goldstino mass. This larger one becomes exactly masslcss 
when the supcrHiggs mechanism is operative as discussed 
in the previous section. This phenomenon can be com- 
pared to the positive complex scalar mass splitting into 
one Higgs boson mass plus a goldstone boson mass when 
the Higgs mechanism is in operation. The correspon- 
dence is what the axino mass to the goldstino is what 
the Higgs boson mass to the goldstone boson mass. 

But there is also a non-gravitational contribution to 
the axino mass. Even if the axion is masslcss, both 
fermion partners are massive as shown in Eq. ([6]) . This is 
a supersymmetric mass where fermionic masses are split 
differently from those of bosons. As an example, consider 
a KSVZ coupling 



W= -f Q SQQ. 



(40) 



There exists a two-loop diagram generating the axino 
mass as shown in Fig. [5] The mass estimate is 



(- 47r 2-)2 TO Sluino - ^ Wgluino- 



(41) 



So, the axino mass contains the contribution 

m~a = ^2 £a m l/2,aH (42) 



a— gaugino 



D. With accidental symmetries 

In [f|, the possiblility of kcV axino was discussed in 
case the superpotential has an accidental symmetry. The 
keV [i| and even eV Q range axino masses are possible 
with some accidental symmetries. The accidental sym- 
metries may forbid the leading order masses of the scales 
m 3/2 and m 1/2 , a - 

In the gauge mediation scenario, to 3 / 2 is negligible and 
the axino-goldstino mixing does not give a significant 
contribution. Then, the loops may give the dominant 
contribution. But the accidental symmetry may forbid 
diagrams of the form of Fig. [5] The superpotential may 
introduce a nonrenormalizable term suppressed by Mp, 
and the expansion parameter is f a /Mp ~ 10~ 7 . Thus, 
the axino mass of Fig. [S]is further suppressed by ~ 10 — 7 
and wc expect 10 GeV ■ 10~ 7 ~ 1 kcV. If it is further 
suppressed, then the estimated axino mass is of order 
10" 3 eV. 

In the gravity mediation scenario, m 3 / 2 is a TeV scale 
and the axino mass depends on the Kahler potential. 
Without the axino-goldstino mixing, it was commented 
in Subscc. Ill C[ and the discussion with some accidental 
symmetries is the same as the above paragraph. 
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III. PARAMETRIZATION WITH 
NONMINIMAL KAHLER FORM 

Chun and Lukas studied the axino mass with the min- 
imal Kahler form (To| . Here we go beyond the minimal 
Kahler form, work with the PQ symmetry realized in the 
NG manner, and include the effects of F terms of the PQ 
charged fields which affect the axino component. 

Equation (|16[) is our definition of axion superfield, and 
the Kahler potential must respect the shift symmetry of 
axion. Therefore, the lowest order terms in the Kahler 
potential with some mixing with SUSY breaking coaxino 
C are 1 



K =l(A + A*) 2 + e(A + A*)(C + C) 
+ CC* + M(A + A*). 



(44) 



The SUSY breaking can be introduced in terms of 
Kim's generalized form [l|| of the Polonyi one [3l| . Here, 
however, we will parametrize the SUSY breaking just by 
an auxilliary holomorphic constant 0, 



9 = 1 + m s tf 2 



(45) 



In view of the discussion of Subsubscc. Ill A 21 if there are 
coaxions then the supcrpotcntial can be taken as 



W(C) 



Wp' 



e 



(46) 



with (W{C)) = M 3 ~ (10 13 GcV) 



A. Local SUSY with one axino 

The most important requirement is that goldstino is 
defined in the vanishing cosmological constant(CC) vac- 
uum, satisfying the U(l) invariance condition. Thus, in 
calculating the axino mass, we satisfy the following three 
conditions: 



(47) 



(i) The vanisihing CC condition, 

where where G = K + Mf, In \W\ 2 . 
(ii) The vacuum stabilization condition, 

G-'C^V.G, + d = 0. (48) 

(Hi) For the U(l) invariance condition, we use 

K = K(A + A*,C, C) 
W = Ge aA/fa W(C). 



(49) 



1 The M term contributes to the saxion derivative in the La- 
grangian as Md 2 s, which has no effect to the axino mass. 



If there are more than one coaxino, we have 

W = W{C) e aA/f " x e aAl/fl x • • • (50) 



The supcrpotcntial in (|49l) preserves the shift sym- 
metry of A since in G = K + \n \W\ 2 , the |IU| 2 part 
is read as |IU| 2 = |IU(C*)| 2 9e a ( A+i *'/^. 

Now, let us calculate the axino mass given by 

m = m V2 [V i G j +'^G i G j ] (51) 
for two classes of (C) = and (C) 7^ 0. 



1. Case for Ga = and G A / 



As an example for Ga = but for G A 7^ 0, we consider 



K = -(A + A*) 2 + CC* + e(A + A*){C + C*) 
W = e aA lf«W{C) 



(52) 



which is an example for Ga = but G A 7^ 0; so 
goldstino has some axino component. From G — K + 
M 2 In \W/M 3 \ 2 , we obtain 



G A = (A + A*) + e(C + C*) + 



a 

Ta 



2 W C 



(53) 



G c =C* +e(A + A*) + M 2 P -^- 



from which the Kahler metric elements and their inverse 
elements are given by 



G AA 
G AA 

G CA 



1, G , 



e, G 



e, G ( 



1 ^AC — fc i u Ci — t > ^CC 

G AC = 

:, G CC 



l-e 2 
c 



l-e 



l-e 2 ' 
1 



(54) 



The Gij elements which are needed for the mass matrix 
elements are given by 

Gaa = 1, Gac = e, 



G 



cc 



Mi 



W C c _ ( W c \ 2 

w v w ) 



(55) 



Since Gq are constants, G^j, — 0. For this reason, the 

Christoffel symbol Y % - k = G^G^f = 0, and hence the 
Kahler geometry is flat. 

The vacuum conditions (i) and (ii) are given by 



G'Hl.G, = 3MJ,, 
G^C.Y.G, + d = 0. 



(56) 



Since we defined the axion superfield A as the supersym- 
metrization of Goldstone boson a, the VEV (A) is zero 
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as discussed in Sec. HU Now, we calculate for two cases 
of (C) = and (C) ^ 0. 

(i) (O=0: 

Three vacuum conditions of Eqs. (|47"1) and (j48|) give 
the following, 

a 



1 



W c 



w 

Wcc 
W 



Ml 
M 2 P - e 2 



(57) 



= 2 



W c /W 

wUw* 



The first of these requires a = 0, and W does not depend 
on A. Thus, in calculating Ga = the In | 2 part 
of G docs not contribute and we cannot constrain the 
supcrpotential from the Ga = condition. 

The axino mass can be obtained from the mass matrix 
in the (A, C) basis as, 



m 3 / 2 



1 

e M 2 P 



Wcc 
W 



( Wg\ 

V w J 



(58) 



The vacuum conditions determine the (22) element to be 
e 2 . 

To make fermion mass term canonical, we have to re- 
define fermion as ip — > Vip where 



V — US 



J_/l -1 
\/2 V 1 1 



AL+e 





(59) 



and the mass matrix m = V T {ViGj + (l/3Mp)GiGj)V 
is given by 



id 



m 3 / 2 



l + e 



i 2 l-e 



(60) 



Note that |e| < 1 is requred to make fermion kinetic term 
positive definite. The eigenvalues of Eq. (f60j) are and 
m 3 /2- This confirms that the vanishing goldstino mass 
comes out right. The coefficient of m.3/2 in the axino 
mass, ^goldstino, is 1. The axino mass independence of 
e shows that goldstino is defined in the direction of Gc 
only. 

(ii) (O^O: 
Here, the vacuum conditions read 



€(C + C*) + —M 2 

Ja 



o. 



l-e 2 



C 



W- 



3Mp, 



Wcc 
W 



Mi 



(W c \ 2 
V W ) 



C* + (Wc/W) 

c + (wz/w*y 



(61) 



The left-hand side of the first equation is Ga, and hence 
SUSY is unbroken in the A direction. But, for C, we 
have 



(C + C*) = - 



a Ml 

e fa 



(62) 



for /„ ~ 10" 7 A/ P . Hence, for the VEV of C staying at 
the Planck scale, a needs to be small at O(10~ 7 ). Now, 
the fermion mass matrix is given by 



m-3/2 



\ 



Ml 



w 



( Wc\ 

\ w ) 

+ M P ^-\ 

3 \M P ^ nlp W ) 



(63) 



where the vacuum conditions are used to simplify the 
(22) element. 

We should redefine fermions such that the kinetic 
terms are canonical, using Eq. (|59j) . The eigenvalues 
are again and TO3/2, confirming the correct goldstino 
mass eigenvalue. The coefficient of m 3 / 2 in the axino 
mass, Cgoldstino, is 1 again. 



2. Case for G A / 



In Case 1, we considered K = \{A + A*) 2 + CC* + 
e(A + A*)(C + C*), which gave G A ^ even though 
Ga = from the vacuum condition. To investigate the 
case of Ga 7^ from the beginning with no mixing with 
SUSY breaking coaxino C, Fig. [H let us consider 



K = f 2 e Cl{A+A ' )/fa + e - c ^ A + A ')/f" 

W = e aA/f "W(C). 
Then, with A = 0, we obtain 



CC* 



(64) 



K A = fa 



Cl (A+A*)/fa _ 



c 2 e 



-C 2 (A+A*)/fa 



fa(ci - C 2 ), 



(65) 



K C = C* 



and 



K AA 



_2ci{A+A*)/f a 



c 2 e ~c 2 (A + A*)/f a 



= (ci + c 2 ), 

Kac = Kcc = 0. 
The Kahler metric is given by 



(66) 



K 



AA 



c 2 e a(A+A*)/f a +( a e - C2 (A+A*)/fa 

(ci + c 2 ) 



K 

K, 



AC 
CC 



K, 



CA 



0, 



(67) 
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and its inverse is given by 
1 



K AA = 



K A A 



K cc = l, K AC = K CA = 0. 



(68) 



Since 

^ AAA 



1 

T, 
i 

Ta 



c 3 eCl (A+A*)/f a 



(4-4), 



3„-c 2 (A+A*)//„ 



(69) 



and other 's vanish, the only nonzero Christoffcl 
symbol is 



AA 



1 ( 4-4 \ 



fa 



(70) 



In fact, this is where the vanishing contribution of the 
superpotential to Gac comes in. Suppose A = C = 
and Ki are negligible as assumed in [32|. Then, d ~ 



M£Wi/W. Let ViGj- = K, 



(W,W 3 /W 2 )} = Xij - Mf,(WiWj/W 2 ). For i = A where 
we are interested in, the second vacuum condition reads 



= G A V A G A + G C V A G C + G A 
w 2 

= G A [X AA -M^] + G C [X 



AC 



M; 



W A W Cl 



p W 2 

~ G A X AA + G c X A c — G A [G A G A + G c Gc]/Mp + G A 

(76) 

where G t ~ MpWi/W is used. Since [G A G A + G C G C ] = 
3Mj> from the first vacuum condition, we have 



G A 



Now, consider the SUSY breaking. The fermion mass 
matrix. Gj, the barometer of SUSY breaking, is 



x ? Ga G ° X 

AAA - l-^X ~ -^jX AC . 



(77) 



G A = / (cx - c 2 



-Mr 



(71) 



Go = C* 



W 



Mi 



And, Gij in the mass matrix is given by 



G AA = (4 + 4) 
G AC = 

-Wcc 



Gcc = Mi 



(W c \ 2 

v w ) 



(72) 



Note that G A c = 0. It is because of the form of the 
superpotential. If the superpotential is merely given by 
a sequestered form, W = W^ a \A) + W^(C), the con- 
tribution of the superpotential to G A c is 



W AC W a W c 



W 2 



w A a) w£ ] 
w 2 



(73) 



so that it does not vanish in general. However, the 
shift symmetry of A in G restricts the A dependence 
of the superpotential to the form W(C) e aA '^ a so that 
(Wac/W) — (Wa Wg /W 2 ) vanishes. This seems not giv- 
ing ma = 2m3/2 of Rcf. [32|]. The reason is the following. 
The vacuum condition G^GiGj = 3Mp reads 

c l > c 2 V Ja ' 



c* 



Ml 



Wc 
W 



(74) 



and the conditions G J ViGj 



3M% 



Gi — read 



(4-4) 
(4 + 4? 

Wcc 



(ci 



. M 2 
c 2 ) + a— 

J a 



M 



W 



/Wcy 
\ w J 



c* 



2, 

-MpWc/W 



(75) 



C + M P W* G IW* 



The factor 2 in front of (Ga/G a ) is the source of 2m 3 / 2 
of HI. For this factor, W A W C /W 2 in V A G C plays a 
crucial role but is vanishing with our potential, and we 
cannot use the first vacuum condition for the factor 2. 
So, the PQ invariance in the superpotential does not re- 



quire rria ~ 2to 3 /2- In fact, Eq. (|76p shows how factor 
2 of Ref. comes about. The factor 2 in the RHS 

of (|77|) is its origin, m = 2m 3 / 2 - Note that it is based 
on the formalism of Eq. (|73p . the result of a sequestered 
superpotential. This does not hold in our case, as ex- 
plained below Eq. ([73"]) . In our case, the superpotential 
of the form e aA ^ a W(C) leads to different relations on 
the mass term and other parameters from those of Rcf. 
[32| . where the Kahler potential does not play a crucial 
role. 

Even though our examples with one coaxino show 
ma — "i3/2j h does not imply that the axino mass should 
be the gravitino mass, as commented later below Eq. 
(|82[) . Moreover, if we put mixing between A and C in 
the above Kahler potential, (1/M%)(A + A*)CC*, the 
axino mass is not the exact gravitino mass. It may be 
larger or smaller than the gravitino mass, but no huge 
enhancement beyond of order 1: 77ia = 0(\) x m 3 / 2 . 

Furthermore, vacuum conditions imply that 



fa< 



\4 



2 (c\+c 2 f/ 2 



21 M P 



(78) 



which would be a criterion for reasonable PQ scale. 

From the vacuum conditions, the fermion mass matrix 
is given by 



rriij 1 
— = ViGj + -G l G ] 

(-(4 + 4) o 



1 / G\ GaGc 
3 y GaGc G 2 c 



(79) 
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where G 



c 



Ce v 



In the G l direction, we confirm 



rriijG 3 = 0, and we have the massless fermion in the 
G % direction, as expected. 

The canonical kinetic terms can be obtained by rescal- 
ing, 

/ I n \ 

(80) 

Then, the physical mass matrix is given by 

m ab ( -1 \ , 1 ( G A G A , G a G c 




-1 

-e 2luJ 



3A/J, V G A G C , GcG 



(81) 

since G A = (1/4 + c\)G A and G c = G c . Then, the 
axino mass Eq. (73) should be modified to 

1 



m 3/2 



-tt(GaG a + GcG c ) 
p (82) 



since G A G A + G C G C = 3A/J,. By the phase redefinition, 
we obtain = 777.3/2- 

However, when c\ = c 2 = c, the situation changes dras- 
tically. The vacuum condition requires Ga = (a/ f a ) = 
so that axino does not break SUSY and is completely 
decoupled from the SUSY breaking field Z = C. So, the 
axino is still massless. It may be a low energy effective 
description o£W= Z 1 {S 1 S 2 - A 2 ) + Z 2 {S 1 S 2 - f 2 ) dis- 
cussed in [l6[ . Even though SUSY is broken in this case, 
the axion superficld docs not appear in the superpotcn- 
tial. Since the axion superfield and the SUSY breaking 
supcrfield are decoupled in both Kahler potential (since 
canonical Kahler is assumed) and superpotential, axino 
is massless since axion is massless. 

Physics of the PQ symmetry breaking can introduce 
another scale of axino mass. Suppose we regard ip as a 
dynamical supcrfiel such that Kahler potential is given 
by, for example, 

1 



K 



e c(A+A*)/f a+e - 



c(A+A')/f a 



(83) 



One may assume that superpotential depends on C and 
if separately, W = e aA/ fa {W\ (ip) + W 2 (C)). In this case, 
we have to consider 3x3 fermion mass matrix in the 
basis (ip, C, A) and find that three eigenvalues are zero, 
0(f a ), and 0(m 3/2 ) + 0(mj /2 /f a ) , respectively. The 
last case corresponds to the axino. When the leading 
order term 0(7773/2) is suppressed from some accidental 
symmetry, the next leading order 0(rriy 2 / f a ) is the ax- 
ino mass scale. 



B. The KSVZ model 

In the KSVZ approach, one introduces the heavy quark 
fields Ql and Q R in the superpotential as j25j . 



KSVZ 



i 3 9 e A l Sa + /qQlQr P e 



A/fa 



(84) 



Model 


Si 


& 


Ql 


Qr 


H d 


H u 


QL 


D C R 


u c R 


KSVZ 


1 


-1 


1 

2 


1 

2 

















DFSZ 


1 


-1 








-1 


-1 


I 


l-i 


1-t 



TABLE I: The PQ charge assignment Q. Ql and Q R denote 
new heavy quark multiplets. 



The PQ symmetry is given near the e point of Fig. [TJ 
with T(Q L ) = -1/2,T(Q R ) = -1/2, and T(X) = 1. 
Near e, there is no ip type field. But near Adw, F(Ql) 
and Qr are not of the ip type, only X is a ip type field, 
and Q obtains the heavy quark mass tuq = fQ(ip(X)). 

It can be rephrased as follows. The heavy quark in- 
teraction with A after integrating out heavy scalars by 
p> = f a , we have the interaction m,QQ[ J Q R e A ^ a . Techni- 
cally, we loose the PQ quantum number information of 
heavy quarks since they do not have a ip type compo- 
nent but only the phase dependence by the original PQ 
charges. These phases can be rotated away by redefining 
the phases of Ql and Qr. This heavy quark interaction 
with A generates the two loop mass presented in Fig. [5] 
at the order of 10 GcV. 



C. The DFSZ model 

In the DFSZ framework, the SU(2) L xU(l)y Higgs 
doublets carry PQ charges and thus the light quarks are 
also charged under U(1)pq [2f3|. The charge assignment 
is shown in Table HI So, the superpotential is written as 



Wdfsz = W PQ + l^-SlH d H u , 

MP 



(85) 



where HdH u = e a pH2H^. Integrating out Si, we have 

Wv ¥S7j = Lie 2A /^ip{H d )p(H u )e- 2A lf* 

+ f u q L e l6 u R e^ s p(H u )e- A /f« (86) 
+ f d q L e is d R e^ e ip(H d )e- A ^. 

For the quarks, they do not contain the ip type fields since 
they do not contribute to V 2 of Eq. (fTi?|) and their phase 
is just a phase parameter 9. This 9 can be removed by 
redefining the phases of quarks, and we obtain 



Wdfsz = A* 



V u V d 



+ (mttLt% + mbbL-bf 



)e 



A/f a 



(87) 



with Adw = 6. Equation (|87|) breaks the PQ symmetry 
and the axion coupling to stop is 



M% In 



(2M6 



fiv u v d m t t L t c R e A/fa + h.c^j 



j^ 2 -/j,v u v d mt\{tLt R }\ cos 



fa 
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whose coefficient is much smaller than the QCD instan- 
ton contribution Aq Cd . So we can neglect the constant 
term for the axion mass even though it breaks the PQ 
symmetry. Also, the constant term does not contribute 
to the axino mass. Including the SUSY breaking auxil- 
liary field 0, we consider 

Wdfsz = I^—tt— + (m 3 Q + m t t L t c R 

N (89) 
+ m b b L b c R + ■■■je A/f ". 

Through the quark mass terms, we obtain the axino mass 
as shown in Fig. [5] with Q replaced by the SM quarks. 
The SM quark mass is at most mt/niQ ~ 10 -9 smaller 
than that of the heavy quark and we obtain the axino 
mass in the range 10 eV. Only the SUSY breaking soft 
mass contribution m 3 can contribute to the axino mass. 



IV. CONCLUSION 

After properly defining the goldstino and axion multi- 
plets, we presented the calculational scheme of the axino 



mass in the most general framework. For only two light 
supcrficlds of goldstino and axino, we obtain £ go idstino 
of Eq. d43j. For G A = where G = K + \n\W\ 2 , 
we obtained ma = m 3 / 2 with the axino- gravitino mix- 
ing parameter e in the Kahler potential. For G a 0, 
we showed that the axino mass depends on the details 
of the Kahler potential. But there is another parameter 
proportional to the gaugino masses, and we can take a 
wide range of the axino mass for cosmological applica- 
tions jlTtlol]. If the gravity mediation is the dominant 
one, the axino mass is probably greater than the grav- 
itino mass, but its decay to gravitino is negligible. Still, 
it softens the cosmological gravitino problem {33j some- 
what as discussed in Ref. H. 
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